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ABSTRACT 

Pulsars are unique astrophysical laboratories because of their clock-like timing 
precision, providing new ways to test general relativity and detect gravitational 
waves. One imp ediment to high-p recision pulsar timing experiments is timing 



noise. Recently iLyne et al.l (120101 ) showed that the timing noise in a number 
of pulsars is due to quasi-periodic fluctuations in the pulsars' spin-down rates 
and that some of the pulsars have associated changes in pulse profile shapes. 
Here we show that a non-radial oscillation model based on asteroseismological 
theory can explain these quasi-periodic fluctuations. Application of this model to 
neutron stars will increase our knowledge of neutron star emission and neutron 
star interiors and may improve pulsar timing precision. 

Subject headings: pulsars: general — stars:neutron — stars:oscillations — white dwarfs 



1. Introduction 

Pulsars have long been known as superb laboratories for fundamental physics due to 
their timing precision and stable pulse profiles. Timing observations of pulsars have re- 
sulted in the best tests of general relativity in the strong-field regime, the most precise 
stellar mass measurements, and the most constraining limits on gravitati onal wave sources 



in the nanoHertz regime through 'pulsar timing array' experiments (e.g. iJenet et al.l 12006 
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Kramer et al.ll2006l ). One significant limitation for all of these studies is timing noise, or 



non-random varia tions in the residuals between measured and model-predicted arrival times. 



Lyne et al.l ( 120101 ) recently showed that the timing noise in 17 pulsars is due to the pulsar's 



spin-down rate fluctuating between different values in a quasi-periodic manner. Furthermore, 
for six of these pulsars, the changes in the spin -down rate directl y correlate to changes in 
the pulse shape. While the pulsars described in iLyne et al.l (120101 ) have long periods, some 
millisecond pulsars, such as PSR B1931-I-27, also exhibit timing noise and may show similar 
but less obvious effects. This behavior offers a possible prescription for removing timing 
noise from data and dramatically increasing timing stability, facilitating projects such as 
the direct detection of gravitational waves through pulsar timing. In this paper, we present 
a model for this behavior based on non-radial oscillations. We outline the model in ^ 
discuss underlyin g assumptions in |J31 describe the implications for the pulsars discussed in 
Lyne et al.l ( 120101 ) in §11 and present conclusions in SJ^l 



A Non-Radial Oscillation Model 



We s how that a non-rad ial oscillation model can explain all the long-period variations de- 
scrib ed in lLyne et al.l (120101 ). In our model, based on astroseismological theory (jPziembowski 
19771 ). non-radial os cillations of high spheric al degree (£) are aligned to the magnetic axis 
of the neutron star (IClemens fc RosenI 120041 ). Between any two adjacent nodal lines, dis- 
placements of material move horizontally along the stell ar surface toward an d away from 
the magnetic pole, as shown in the schematic of Figure [H iJessner et al.l (120011 ) showed that 
thermal emission of electrons dominates the rate at which charged particles flow from the 
surface into the magnetosphere. Non-radial oscillations, and their associated velocities, will 
vary the amount of thermal emission by creating regions of local heating due to compression 
between the nodal lines, thus increas ing the flow of cha rged particles into the magnetosphere 
and modulating the radio emission (IStrohmayerlll992l ). 



The situation is analogous to behavior seen in white dwarf stars. Some white dwarfs 
are g'-mode, multi-periodic pulsators that tend to excite a select number of pulsation modes, 
as shown in Figure |2l In adiabatic oscillations, local heating of the stellar surface coincides 
with pulsation maximum during the compression phase of the pulsation. A quarter cycle 
later, material will be moving away at a maximum velocity. The displacements responsible 
for surface heating and flux chang es in oscillating whit e dwarfs are primarily non-radial 
because of their high surface gravity (iRobinson et al.lll982l ) and ha ye measured radial velocity 



variations that can b e interpreted as the horizontal velocities ( Ivan Kerkwijk et al.l 12000 



Thompson et al.ll2003l ). In neutron stars, we expect that the amplitude of the displacements 
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Fig. 1. — Schematic: A map of the non-radial oscillations on a pulsar. The circle around 
the magnetic pole in the enlarged view denotes the boundary of the emitting region. Dark 
regions indicate negative displacements and light regions positive ones. After a half cycle 
of the oscillations, the dark regions would be light and vice versa, but the nodal lines (as 
defined by £) separating them would remain unchanged, except for rotation of the whole 
pattern about the rotation axis of the star. The magnetic cap is crossed by four representative 
sightlines, where /3 is the angle between our sightline and the magnetic cap. The listed values 
for /3 are equally spaced. Panels: For each sightline, the columns show the corresponding 
average profile (i.e. intensity vs. pulse phase) for various values of Vosc and i. 
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foUows surface thermal variations caused by non-radial oscillations of the neutron star. As 
with the white dwarf stars, an increase in the velocities can result in increased thermal 
emission from the star. 

In our non-radial oscillation model, each frequency of oscillation has its own peri- 
odic velocity (or displacement) amplitude. In the pulsars we h ave studied to date, the 



amplitude of the displac ements dictates the average pulse shape (IRosen fc Clemens! 12008 



Rosen fc DemorestI 120101). but as the ampl itude of the velocities grows, the average pulse 
shape can change (jClemens fc Rosenll2008l ). In these pulsars, through fitting a non-radial 
oscillation model to their drifting subpulses, we have also found I to be large compared 



to that seen in white dwarf stars, consistent with thei r smaller size. For 
and B09 43+10, ^ values range f rom 126 < d, < 133 (IRosen fc Demorest 



pulsars B0809+74 



2010) and 385 < 

< 875 (iRosen &: Clemensll2008l ). respectively. In white dwarf stars, modes with high £ fade 



from view b ecause of geometric cancellation of the stellar su rface, leaving only modes of low 
£ <4) jYeates et all I2OO5I : [Thompson et al.ll20oi . l2008h . 



By fitting a non-radial oscillation model to the drifting subpulses in PSR s B0943-I-10 and 
B0809-I-74, we measured the pu lsation period to be on the order of 30—50 ms (IRosen fc Clemens 



2008; Rosen fc DemorestI 



20101). These values for the pulsation period are consistent with 



core g'-modes (iReisenegger fc Goldreichlll992l ). However, core q -modes require la rge exci- 
tation energies, and therefore surface g'-modes are a possibility (IStrohmaverl 119931) as they 



have lower energies and larger surface amplitudes than the core modes (IMcDermott et al. 



19881 ). The peri od predictions for core ff-mode s in neutron stars range from a minimura 
value of 10 ms (IReisenegger fc GoldreichI 119921 ) to 2—88 seconds (IMcDermott et al.lll988l ) 
depending on the model for the structure and co mposition of the stellar interior. Surface 
g'-modes have periods in the 40—400 ms range (IMcDermott et al.lll988l ). While the pul- 
sation periods calculated for most models a re calculated assuming low £, the period of the 
mode decreases for higher spherical degree (IMcDermott et al.lll988l ). It is unlikely that the 
pulsation modes we see in PSRs B0943+10 and B0809-I-74 are p-niodes as these modes have 
periods on the order of tenths of milliseconds ( IMcDermott et al.lll988l ) and have overtones 
with shorter periods, whereas the overtones of g-modes have longer periods. 

The two oscillation driving mechanisms for white dwarf stars are the k — 7 mecha- 
nism and convective driving. For pulsating DOVs stars like PG1159— 035, characterized 
by an atmosphere that lacks hydrogen but shows helium, carbon, and oxygen absorption 
lines, the ft: — 7 mechanism drives the pulsations as the opacity varies steeply with pres- 
sure Jstarrfield et al.lll983l : ICorsico et allbood : [Co^l2003h . In DAV and DBV stars, which 
have atmos pheres domina t ed by hydrogen or helium, respectively, conve ctive motions drive 
pulsations (iBrickhilll Il99ll : IColdreich fc Wul Il999l : IWu fc GoldrdchI Il999[ ) . 
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Fig. 2. — The flux variatio ns of GD358, a helium atmosphere pulsating white dwarf star, 
observed on different nights (jProvencal et al.ll2009t ). For two nights in 1996, the light curve of 
GD358 drastically changed. Unlike during all other observations of this star, all the pulsation 
energy was pumped into the shortest period modes. Reproduced by permission of the AAS 
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Driving mechan isms in for neutron star oscillations are different and less well understood. 



Jessner et al.l (l200l[ ) show that thermal and field emission from the neutron star surface can 



accelerate electrons along open field lines with the formation of a vacuum gap. Therefore, 
oscillations should be related to thermal variations at the neutron star surface, as would 
be produced by g'-modes. The high spherical degree that we observe in the data (and is 
predicted by our model) suggests that the pulsation driving energy is concentrated in a 
small area on the surface of th e star. Otherwise it would average away over multiple surface 
zones (jClemens fc Rosen! 120041 ). In one possible scenario, the source of the driving energy is 
at the magnetic cap and the particle emission exerts a torque on the open field lines, which 
is coupled to the rest of the star via magnetic and mechanical dissipation. If the magnetic 
field is coupled (even weakly) to the surface and it displaces material laterally, the torque 
on the open field lines results in a displacement of material on the polar cap. The heating 
that results can increase particle era ission, which then increases the torque, and drives the 
oscillations ( IClemens fc RosenI 120041 ). 



^Other possibilities for dri ving mechanisms iiiyolvin g direct shaking of field li nes (IBoriakofl 

19761 ). thermal stratification (IMcDermott et al.lll983l). chemi cal stratification (jFinru 119871 ). 
and stratification in the core (IReisenegger fc Goldreichlll992l ). The gravitational radiation 
dampi ng timescale for q-modes exceeds the age of the Univer se and is therefore a negligible 
effect ( IMcDermott et al.lll983l : IReisenegger fc Goldreichlll992[ ). The pulsations we observe in 
PSRs B0943+10 and B0809+74 appear to persist on hour to month (if not longer) timescales 
and any driving mechanism must be in agreement with these observations. Other neutron 
star oscillations have been detected in the form of quasi-periodic oscillations (QPOs) from 
soft gamma-ray repeaters. These have been interpreted as torsional shear modes; th e fre- 
quencies of QPOs are too high to originate from g'-modes (jWatts fc Strohmayerl 120071 ). 



3. Standard Model Assumptions 



In the standard pulsar model f lGoldreich fc Julian Ill969[ ). the magnetosphere is consid- 
ered to be 'force-free' in that the electromagnetic energy density dominates all other dissipa- 



tiye forces and that radiative emission has a negligible effect on the spin-down rate (lArons 



20071 ). Instead, the configuration of the magnetosphere, based on the polar cap cascade zone 
and the current densit y within the op e n field lines, and its change over time di ctates the 
pulsar spin-down rate ( iTimokhinI 120061 : lArons 1 120071 : iTimokhiru 120101 ). Recently, iTimokhin 
( I2OIOI ) speculates that the magnetospheres of some pulsars can have several quasi-stable 
states with different configurations and that the switching of th e magnetosphere between 
these states can result in the observed mode changes and nulls. iTimokhinI ( I2OIOI ) further 
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notes that each state would have a different spin-down rate. 



Timokhinl ( 120071 ) shows that oscillation modes of high spherical degree and non-vanishing 



veloc ities at the surface can alter the local Goldreich- Julian (GJ) charge density (iGoldreich &: Julian 



19691 ). The GJ charge density, while not a contributing factor to the spin-down rate in the 
force-free model, is then a combination of the charge density due to rotation and the charge 
density due to oscillations. In this paper, we assume that, like white dwarf stars, neutron 
star oscillations can change with both amplitude and spherical degree, therefore affecting the 
total charge density and the magnetospheric configuration. If the charge density is at least 
partially responsible for the energy loss of the pulsar and/or the change in pulsation mode 
changes the magnetospheric configuration, we expect to observe a change in spin-down rate. 



Theoretical Predictions 



Under some models we expect for the spin-down rate and the pulse shape to c hange with 



both s pherical degree and the amplitude of the velocities of the oscillations, Vosc- iTimokhin 
(120071 ) calculates the GJ charge density for neutron star oscillations assuming the standard 
pulsar model of a rotating magnetized conducting sphere surrounded by plasma, and ex- 
plores the excitation of oscillations by neutron star glitches. He shows that oscillation modes 
of high spherical degree and non-vanishing velocities at the surface can alter the local GJ 
charge density (pcj) where pcj ~ i{Vosc/c){B/4:TTRi\is)- The charge density will affect the 
accelerating electric field. For the oscillations to be st rong enou g h to i nfiuence the particle 
distribution, resulting in effects that we can observe, iTimokhinI ( 120071 ) suggests the spher- 
ical degree is a pproximately several hundred, which is in agreern ent with our independent 
measurements (IRosen fc Clemens! l2008l : iRosen fc DemorestI |2010| ) . Because more particles 
leave the stellar surface when the charge density is greater, we expect the spin-down rate to 
increase with either: 1) increasing i, 2) increasing Vosc, or 3) increasing i and V^.^r;. Assuming 



that the current density in the magnetosphere is given by Equation [T] (ILorimer fc Kramer 
2004 ). the charge density is directly related to the spin-down rate by: 



Pgj 



VLB 
2nce 



7 X 10^°cm"^ 



-1/2 



P 



1/2 



10-15 



where Q is the angular rotation frequency, B is the magnetic field, P is the spin period of the 
pulsar, and P is the spin-down rate. Note that this assumption is not in agreement with the 
simple force-free model. We calculate the fractional change in charge density using the change 
in spin-down rate and then use this to calculate P2/P1 = ^2Ksc,2/^iKsc,i- For the pulsar with 
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the largest change in spin-down rate, PSR B1931+24, the value of iVosc will change by 20%; 
for PSR B1929+10, which has the smallest change in spin-down rate, iVosc will change by 
0.15%. To make these estimates, we assume that the contribution to the charge density due 
to oscillations is much greater than the contribution due to rotation. Figure [T] illustrates 
how changes in these parameters affec t the pulse shapes and Table [1] shows our best guesses 
for the classifications of the pulsars in iLyne et al.l ( l2010l ) into the different scenarios. 



In the first scenario where Vosc changes but i remains constant, pulsars fall into two 
categories: small and large velocity oscillations. For small velocity oscillations, the change in 
the spin-down rate is less than some value, APcrit, and is accompanied by little o r no change 



i n the pulse shape. With the small sample of pulsars that show this behavior in lLyne et al. 



(I2OIOI ). we cannot determine the precise demarcation point between large and small changes 
in the spin-down rate, but estimate it to be APcrit ~10% by inspection of those pulsars 
which show pulsar shape changes. The pulsars in the first column of Table [1] fall into this 
category. Many stable pulsating white d warf stars have subtle changes in amplitude over 
time, thereby changing Vosc but not i (jProvencal et al.l l2009l : IVuilld I2OOOI : [Handler et al. 



20081 ). When the amplitude of a single mode changes, both its flux and velocity amplitudes 
change. If amplitude changes occur on a neutron star, the change in Vosc will also change 
the spin-down rate. A periodic change in amplitude will result in a periodic change in the 
spin-down rate. Because the spherical harmonic is responsible for the average pulse shape 
and i is constant, the average pulse shows little or no change with varying spin-down rate. 

Pulsars with large velocity oscillations have changes in the spin-down rates greater than 
APcrit with changes in the pulse shape, such as B2035+36, are listed in the second column of 
Figure [H In the small sample of pulsars we have quant itatively fit, the average pulse shape 
is dominated by the amplitude of the displacements (IClemens fc RosenI |2008| ). However, 
when the relative amplitudes of the velocities and displacements becomes dominated by the 
velocities, the pulse shape c an change and we may observe a corresponding change in the 
polarization behavior as well (jClemens &: Roserul2008l ). PSR B2035-I-36 shows a sudden large 
change in spin-down rate. Analogously, the white dwarf star GD358 suddenly pumped all of 
its energy into the shortest period modes, increasing the overall amplitu de and changing the 



pulsation characterist ics on a timescale of hours as shown in Figure [2] (IKepler et al.l 12003 



Provencal et al.l |2009| ) . Such exchanges i n the energy between pulsations on white dwarfs is 



not a lways so dramatic but are common (IKleinman et al.lll998 



2OOII). 



Vuillell2000l : IWu fc Goldreich 



Pulsars that have changes in the spin-down rate less than APcrit accompanied by 
changes in the pulse shape, like those listed in column three in Table 1, fall into the second 
scenario, where £ varies but Vosc remains constant, as illustrated in the third column of 
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will modulate both the spin-down rate, the pulsation period 
White dwarf s t ars are known to variably e x cite pulsations 



Figure [TJ Periodic changes in 
and the average pulse profile 
of varying from 1 < i < 4: (iKepler et al.l 120031 : [Thompson et al.l 120081 : iClemens et al. 
20001 ) ■ When oscillations of different i are excited on a pulsar, the relationship between 
the nodal structure to the pulse width is more complicated than a simple proportion due to 
the geometry of the star and our viewing angle relative to the magnetic pole. While it is 
easy to assume that as the area between nodal lines gets larger the pulse width increases; 
for certain geometries, the pulse width can actually decrease, as in PSR J2043-I-2740 and as 
illustrated in Figure [TJ Radio polarization data would allow us to determine the geometry 
and the amount by which i is changing. 

The third scenario, where both Vqsc and i change, can also be possible through the 
transfer of energy between modes. This is shown in the fourth column of Figure [T] and 
Table [H If both Vosc and i vary, we expect larger changes in the spin-down rate and pulse 
width, like that seen in PSR B193H-24. If Vosc and i are changing at different rates, 
we also expect to see more complicated behavior including beat frequencies which would 
result in complicated patterns in the spin-down rate and pulse shape. Extrapolating the 
more complex behavior of unstable white dwarf pulsators to neutron stars, we would expect 
to see multiple, possibly simultaneous pulsation modes. The changes in the light curves in 
white dwarf stars are not necessary gradual or sinusoidal and therefore we do not necessarily 
expect gradual changes in the spin-down and pulse shape in neutron stars. The velocity 
and spherical degree could change depending on which mode is dominating the pulsation 
spectrum. 

All the scenarios above can explain how the spin-down rate can change periodically with 
or without accompanying changes in the pulse profile. It may be that some pulsars do not 
exhibit any changes in the spin-down rate with time but still show pulsations. These objects 
would have stable pulsation modes. Stable white dwarf pulsators can have seve ral pulsation 
modes present simultaneo usly and tend to be dominated by a few pulsations (lYeates et al. 
20051 : iMuUally et al.ll2008l ). We speculate that neutron stars which do not exhibit changes in 
spin-down rate or period show no changes in Vosc or i, but may have several pulsations modes 
present. PSRs B0943-I-10 and B0809-I-74 have both been quantitati vely fit to an oscillation 



model but do not exhibit any known changes in the spin-down rate (iRosen fc ClemensI 12008 



Rosen fc Demorestll2010l ). Just as only a percentage of white dwarf stars exhibit pulsations, 
we expect only a fraction of neutron stars to show oscillations. 
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Conclusion 



In summary, we have shown that the fluctuations in spin-down rate described in lLyne et al. 



( 120101 ) can be explained by a non-radial oscillation model. Our model is not complete, as 
little is understood about pulsar magnetospheres and there is no clear mechanism for driving 
the oscillations. However, our oscill ation model does not requ ire invoking an external influ- 



ence to explain the periodicity, as in lCordes fc Shannon! ( l2008l ). and is based on observations 



of other well studied stars. This will lead to better understanding of pulsar emission physics 
and may enable optimal removal of effects due to non-radial oscillations in pulsar timing 
data, increasing the usefulness of pulsars as fundamental physics laboratories. 
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and the National Science Foundation. The National Radio Astronomy Observatory is a facil- 
ity of the National Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. Also we thank Duncan Lorimer and Xavier Seimens for helpful comments 
on the manuscript. 
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Scenario 1 



Scenario 1 



Scenario 2 



S 



cenano 6 



AP< APcrit AP> APcrit AP< APcrit AP> APcrit 

no pulse shape changes small pulse shape changes pulse shape changes large pulse shape changes 

large A^,,, no Ai no A^.c, Ai A^sc, Ai 



small AVosc, no A£ 



B1642-03 
B1903+07 
B1839+09 
B2148+63 
B1818-04 
B0950+08 
B1714-34 
B1907+00 
B1826-17 
B0919+06 
B1929+20 



B2035+36 



J2043+2740 
B1822-09 
B1540-06 
B1828-11 
B0740-28 



B1931+24 



Table 1: List of pulsars which show changes in the spin-down rate with possible classifications 
with respect to non-radial oscillation properties. 



